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Abstract 

A planar straight-line graph which causes the non-termination Ruppert's algorithm for a minimum 
angle threshold a ^ 29.5° is given. The minimum input angle of this example is about 74.5° meaning 
that failure is not due to small input angles. Additionally, a similar non-acute input is given for which 
Chew's second algorithm does not terminate for a minimum angle threshold a ^ 30.7°. 

For a non- acute planar straight-line graph, Ruppert's algorithm produces a conforming Delaunay trian- 
gulation composed of triangles containing no angles less than a. Ruppert proved the algorithm terminates for 
all a ^ 20.7° |3] and a minor addition to the analysis extends the results to input with all angles larger than 
60°. In practice, the constraint a ^ 20.7° has been seen to be overly conservative. Ruppert observed that 
the minimum angle reaches 30° during typical runs of the algorithm. Further experimentation by Shewchuk 
[4] suggested that even higher values are admissible: "In practice, the algorithm generally halts with an 
angle constraint of 33.8°, but often fails [at] 33.9°." In this note, we demonstrate an input (the upcoming 
Example 2) for which Ruppert's algorithm does not terminate for some minimum angle parameter a less 
than 30°. We begin by revisiting the best known example which causes non-termination for any a > 30°. 

Pav Example Steven Pav gave an example demonstrating that Ruppert's algorithm can fail to terminate 
for any a > 30° [2]. This example, depicted in Figure [1] involves two adjacent segments with lengths 1 and 
-\/2 such that they form a triangle with a 30° angle. Pav observed that the circumcenter of this triangle 
lies on the boundary of the diametral ball of the longer segment. This causes the longer segment to split. 



*This note is an extension of "On the Termination of Ruppert's Algorithm" which appeared in the Research Notes of the 19th 
International Roundtable, 2010. While generated after the submission deadline, these results were presented at the conference 
with the original research note. 
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Figure 1: The Pav Example demonstrates that Ruppert's algorithm may fail for any a > 30° [2]. 
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Figure 2: Example 1 demonstrating non-termination of Ruppert's algorithm (and Chew's second algorithm) 
for a > 30.7°. 

yielding a similar configuration in which the adjacent segments are smaller by a factor of l/\/2 and then 
repeats indefinitely. Note: the point on the boundary of the diametral ball does not technically encroach 
the segment but small perturbations of this configuration yield non-termination of the algorithm for any 
a > 30°. 

Improvements to the Pav Example face a key obstacle: each split must reduce the length of the shortest 
segment in the mesh by a fact of 1/ which is the largest possible reduction according to the theory. To 
eliminate this problem, we consider examples containing more adjacent input segments which allows for a 
smaller reduction in segment length at each midpoint insertion. However, throughout these examples we 
maintain our requirement that input angles are larger than 60°. 

Example 1 Consider a non-acute input containing four adjacent segments of lengths 2, 2'^/'*, 2^/^ and 2^/^ 
as in Fig. [5J The endpoints of the longest and shortest segments form a Delaunay triangle with smallest 
angle arctan2~'^/'* « 30.7°. If a > arctan2~^/^, the circumcenter of this triangle encroaches upon the longer 
segment causing the midpoint of the longest segment to be inserted. Now the adjacent segments have lengths 
1, 2'^/'*, 2^/^ and 2^/"* and the ratio of the shortest and longest segment is still 2^^'^. Again this gives a poor 
quality triangle and the midpoint of the longest segment is inserted. This cycle repeats indefinitely. 

While the Pav Example provides a sharper limit on the performance of Ruppert's algorithm, Example 
1 has importance of its own. Chew's second algorithm [T] (widely used as the default mesh generator in 
Triangle [5j) terminates on the Pav Example while it fails to terminate on Example 1. 

Following the approach of Example 1, input with three and five adjacent segments can be constructed 
seeking non-termination; see Figure [H In the three segment version the skinny triangle circumcenter does 
not encroach the longest segment and no cascading encroachment occurs. The five segment variant does lead 
to non-termination but only for an angle threshold larger than about 33° and thus does not improve upon 
the previous results. 

Example 2 A better example can be generated by combining ideas in the Pav Example and Example 1. 
First consider a four segment input that contains two adjacent copies of Example 1 as shown in FigureHJ^left). 
Thanks to a fortuitous angle of about 29° (at the top of the figure) non-termination can occur for any a > 30° 
by spiraling around all four input segments rather than just alternating between two segments as in the Pav 
Example. 

Finally, improvement is made by perturbing this example to balance the 30° and 29° angles. As shown 
in Figure [Upright) the length of the upper segment (2a) as well as the angle between two segments {9) are 
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Figure 3: Example input with three and five adjacent end segments. Left: The three segment version 
does not lead to indefinite refinement. Right: the five segment variant leads to a weaker restriction on the 
minimum angle threshold. 




Figure 4: Left: a combination of the input in the Pav Example and Example 1. In this configuration, 
Ruppcrt's algorithm still fails to terminate for a > 30°. Right: a generalization of the example to the left. 
6 and a can be varied to minimize the maximum of ai and Q2 while requiring that an infinite encroachment 
sequence occurs. 
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Figure 5: Best known example for which Ruppert's algorithm fails to terminate without small input angles. 

allowed to vary while still requiring that the desired encroachment occurs. Our goal is to minimize the larger 
of the small angles in the skinny triangles (ai and a2 in Figure[4jright)). The minimum value of max(Q!i, 02) 
occurs at a solution of the following four equations. 

sin 6 = cos 9 + a/V2 
cos 9 = 2a — cos ai \/ \d?- + 1 — 4a cos 9 
sin 6* = tanQ;2(cos6' + ^Pija) 
a\ — Q!2 

Searching near the original values 9 = 75°, a = 1, ai = 29°, and a2 = 30° gives a numerical solution: 

61 = 74.51°, ai,2 = 29.51°, and a = 0.985. 

Thus for this configuration, depicted in Figure[Sl Ruppert's algorithm fails to terminate for any a ^ 29.51°. 
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